UC'RL-^OOII 


UNIVERSITY  of  CALIFORNIA 


LIVERMORE 


LOVELACE  FOTJNOAflO^-^ 

D  (.)  C  ■  U  i  t  -  '  i  ^ 


28872 

SEP  5  iqK7 


TID-4500,  UC-35, 
Nuclear  Explosions  — 
Peaceful  Applications 


Xadi'virvexxo  e  XKia.d±cb'fcl.oxi.  XiCi.'b03?Ei>'fc03?3r 

UNIVERSITY  OF  CALIFORNIA 
LIVERMORE 


UCRL-50211 

NECESSARY  CONDITIONS  FOR  SIMILAR  SOLUTIONS 
OF  PROBLEMS  OF  TURBULENT -GAS  DYNAMICS 

Barbara  K.  Crowley 
January  31,  1967 


Reproduced  From 
Best  Available  Copy 


DISTRIBUTION  STATEMENT  A 

Approved  for  Public  Release 
Distribution  Unlimited 


Printed  in  USA.  Available  from  the  Clearinghouse  for  Federal 
Scientific  and  Technical  Information,  National  Bureau  of  Standards, 
U.  S.  Department  of  Commerce,  Springfield,  Virginia  22151 
Price:  Printed  Copy  $3.00;  Microfiche  $0.65. 


CONTENTS 


Page  No. 

ABSTRACT .  ^ 

INTRODUCTION .  ^ 

SECTION  I:  DERIVATION  OF  SIMILAR  CONDITIONS  FOR  THE  PUFL 

EQUATIONS .  ^ 

Continuity .  ^ 

Momentum  2 

Energy  3 

Derivation  of  the  Similar  Conditions  3 

SECTION  II:  AUXILIARY  EQUATIONS  USED  BY  PUFL  ....  7 

Equation  of  State .  7 

Wall  Shear  Stress .  8 

Ablation  Equations .  8 

SECTION  III:  CALCULATIONAL  VERIFICATION  OF  THE 

SIMILAR  CONDITIONS  .  H 

SUMMARY . 15 

ACKNOWLEDGMENTS . 18 

REFERENCES . .  •  18 

APPENDIX  A  19 

APPENDIX  B  22 

Momentum  22 

Energy  23 


-1- 


NECESSARY  CONDITIONS  FOR  SIMILAR  SOLUTIONS 
OF  PROBLEMS  OF  TURBULENT -GAS  DYNAMICS 


ABSTRACT 


The  quasi -one -dimensional  conservation  equations  with  friction,  heat,  and  mass 
sources  are  considered.  Necessary  conditions  for  similar  solutions  of  these  equations 
are  mathematically  derived,  and  these  conditions  are  verified  by  sample  computer 
calculations. 


INTRODUCTION 


PUFL  is  a  quasi -one  “dimensional,  ’’ almost -Lagrangian"  computer  code  for 
solving  the  partial  differential  equations  of  mass,  momentum,  and  energy  conservation 
for  gas  flows  in  pipes. ^  The  differential  equations  are  cast  in  the  Lagrangian  form, 
but  mass  sources  are  allowed.  Also,  frictional  forces  and  energy  sources  or  sinks 
are  considered. 

In  addition  to  providing  a  better  understanding  of  the  equations,  there  are  two 
major  reasons  for  seeking  similar  solutions  of  the  PUFL  equations.  First,  the  number 
of  calculations  that  need  to  be  considered  could  be  reduced;  and,  secondly,  knowledge 
of  the  mathematical  conditions  for  similarity  could  provide  a  basis  for  designing  scaled 
experiments. 

In  Section  I  of  this  report,  the  general  conservation  equations  for  pipe  flows  with 
ablation  are  given  and  briefly  discussed.  The  conditions  necessary  for  similar 
solutions  are  then  derived.  The  auxiliary  equations  used  by  the  PUFL  calculation  to 
simulate  ablation  are  presented  in  Section  II.  The  necessary  conditions  that  the 
similarity  conditions  from  Section  I  impose  on  the  auxiliary  equations  and  on  the  physical 
parameters  are  discussed.  Section  III  presents  some  sample  PUFL  calculations  which 
indicate  that  by  using  the  derived  conditions,  similar  solutions  are  obtained. 
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SECTION  I: 

DERIVATION  OF  SIMILAR  CONDITIONS  FOR  THE  PUFL  EQUATIONS 


The  conservation  equations  that  are  the  basis  of  the  PUFL  model  are  presented 
below.  The  detailed  derivation  of  these  equations  may  be  found  in  Ref.  1. 


Continuity 


The  continuity  equation  is 


Dm  _  D 


Dt 


Dt 


(pV)  =  Sm, 


where  rn  is  the  mass  flux  (mass/ area-time)  entering  a  volume  homogeneously  along 
the  side  walls  of  the  pipe  that  has  a  surface  area  S.  The  entering  mass  is  assumed  to 
mix  instantaneously  with  the  material  already  present. 

For  extensive  (mass -dependent)  variables  such  as  the  volume,  the  substantial 
derivative,  D/Dt,  is  simply  d/dt.  Hence,  after  rearranging,  the  above  equation 
becomes 


Dp  _  0  .  dV  ,  Sm 

Dt  ”  "V  dt  V  • 


Momentum 


The  equation  for  the  conservation  of  momentum  is 
(mu)  =  SF, 

where  SF  is  the  vector  sum  of  the  forces  acting  on  an  element.  Since  the  mass  of  a 
particle  with  mass  sources  varys  with  respect  to  time,  the  mass  is  inside  the  derivative 
operator. 

The  one  “dimensional  momentum  equation  used  by  PUFL  is 
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The  first  term  represents  the  adjustment  of  momentum  in  the  gas  due  to  the  mass  flux 
(m).  The  velocity  of  entering  material  is  u^  in  an  x  direction  (u^  may  be  positive  or 
negative).  The  accelerating  force  due  to  pressure  differentials  is  denoted  by 
The  term  T  S  is  the  retarding  frictional  force  exerted  on  the  gas  flow  by  a  pipe  wall 
having  a  sur^ce  area  S.  Viscous  forces  within  the  gas  are  assumed  to  be  zero.  Also, 
body  forces  are  assumed  to  be  negligible. 


Energy 


The  conservation  of  total  energy  -  kinetic  plus  internal  -  is  considered  by  PUFL. 
By  using  the  momentum  equation  for  velocities,  the  total  energy  equation  reduces  to  an 
equation  for  specific  internal  energy: 


De  _  L 

Dt  m 


2 


2 

~  +  (e 
'  w 


dV 

-PdT 


+  T  S 
w 


+  H) 


The  first  term  represents  the  modification  of  specific  internal  energy  within  a  gas  by 
the  total  energy  of  the  entering  mass.  The  expression  p^  is  the  usual  work  term, 
and  T  S  ju|  expresses  the  conversion  of  the  translational  kinetic  energy  of  a  gas  into  in¬ 
ternal  energy  of  heat  due  to  the  shear  stress  (T,^)  at  the  wall.  The  rate  of  energy 
change  within  a  gas  due  to  sources  or  sinks  is  denoted  by  H .  This  term  be  related 
to  the  mass  entering  the  gas  by  using  an  additional  ablation  equation. 


Derivation  of  the  Similar  Conditions 

In  order  to  obtain  relationships  necessary  for  similar  solutions,  the  conservation 
equations  just  described  are  now  examined.  The  approach  used  in  the  following 
discussion  is  equivalent  to  the  familiar  approach  frequently  used  in  dimensional  analysis 
to  nondimensionalize  flows  by  the  use  of  characteristic  or  reference  parameters.  See 
Appendix  A  for  verification  that  this  approach  is  equivalent. 

2  , 

A  dimensional -analysis  approach,  using  the  Vaschy -Buckingham  vr  theorem,  also 
yields  the  same  results  as  those  derived  here.  Although  the  approach  used  here  is 
more  laborious  than  the  tt  method,  it  tends  to  impart  a  feeling  of  being  more  straight¬ 
forward  when  applied  to  a  number  of  relatively  unfamiliar  equations  and  parameters. 

Consider  two  gas  flows,  both  of  which  can  be  simulated  by  the  PUFL  equations. 
Let  one  flow  be  described  by  a  set  of  unprimed  variables,  the  other  by  a  set  of  primed 
variables.  The  flows  are  examined  at  two  instances  when  comparison  of  the  two  flows 
is  to  start.  At  these  instances,  ratios  are  obtained  for  all  of  the  pertinent  variables 
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and  parameters  in  the  PUPL  equations.  These  dimensionless  ratios,  which  relate  one 
flow  to  the  other,  are  denoted  by  subscripts  and  are  defined  as  follows: 


Axial  distance: 

^0=  F 

Radial  dimensional: 

Density: 

''o  ‘fr 

Pressure: 

Po  =  ^ 

Velocity: 

Specific  internal  energy: 

®0 

Volume: 

Surface  area: 

s  =-^ 

°0  S' 

Mass: 

“o'  iF 

Wall  stress: 

X 

T  •  =  W 

Rate  of  energy 
deposition: 

0  JJ, 

Dimensionless  coefficient 
of  friction: 

CfO  ' 

^  Cf 

Flux  of  mass  entering 
the  flow: 

m 

Specific  internal  energy  of 
entering  mass: 

=  fw 

^wO  e' 
w 

Time  intervals  after  the 
instant  at  which 
comparison  starts: 

Specific  energy  of 
vaporization: 

E 

TT  -  V 

vO 

Dimensionless 
coefficient  of  heat 
transfer: 

C 

^HO 

Turbulent-transpiration 

coefficient: 

„  =JL 
r,' 

Ratio  of  specific  heats: 

7 

U  y 

Velocity  of  entering 
mass: 

il 

o 

Derivatives  for  an  arbitrary  dependent  variable  (cj))  in  the  two  flows  are  related 
by 

9c|)  _  9(()  ^  _  8x*  ^  8<{)  _  1  ,  8<}) 

8x  8x’  8x  8x  8x*  8x* 

and  by 

^  =  M  ■  Dt»  ^  Dt^  D({>  ^  1  ,  D<l> 

Dt  DV  Dt  ”  Dt  Dt»  t^  DV  ' 

The  above  relationships  are  now  used  in  the  conservation  equations.  For  con¬ 
tinuity,  this  means  that 


Dp  „  P  d V  .  Sm 
Dt  ■  V'dF^  —  ' 

^0  Uel  =  ^0^  VO  dV'  ^  ^0 

to  to  -  dt'  VoV 
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and 

Do'  -  o'  I  S'm' 

Dt'  'V'  '  dt'  PqVq  ■  V 

S^Mn^n  .  . 

Hence,  if  ^  =  1,  the  continuity  equations  for  two  physical  systems  are  similar, 

Pq^O 

One  interpretation  of  S^M^tQ  =  PqVq  is  as  follows.  Suppose  that  two  flows  are 
observed  over  time  intervals  t  and  t'  respectively.  If  the  ratio  of  the  change  in  mass  to 
the  original  (or  final)  mass  is  the  same  for  both  flows,  then  SQMQt^  =  Pq^O’ 
flows  satisfy  the  same  continxiity  equation  during  the  respective  time  intervals. 

If  one  is  familiar  with  the  technique  just  used  on  the  continuity  equation,  the 
conditions  for  similarity  may  be  written  directly  by  merely  inspecting  the  equations. 
Appendix  B  gives  the  algebraic  details  for  the  derivation  of  similar  conditions  for  the 
momentum  and  energy  equations.  The  conditions  required  for  similarity  by  the  con¬ 
servation  equations  are  given  in  the  following  array: 


Continmty 


Momentum 

(2) 

(3) 

(4) 

(5) 

Energy 

(6) 

(7) 


U, 


wO 


=  1 


Mq 

^p^p^p 

“^wP^P^P 

UqMo 


=  1 


=  1 


u. 


=  1 


=  1 


(8) 


u 


wO 


=  1 
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(9) 

^wO  _  ^ 

®0 

(10) 

^0^0*0 

Mq 

(11) 

Po^o  . 

(12) 

^0^0 

(13) 

If  the  auxiliary  equation  Mq  =  PqVq  is  used,  conditions  (1),  (3),  and  (10)  are 
identical.  With  conditions  (2)  and  (6),  conditions  (7)  and  (8)  are  superfluous.  Also, 
using  condition  (6),  conditions  (5)  and  (12)  are  identical.  If  conditions  (11)  and  (6)  are 
used  in  condition  (4),  the  latter  reduces  to  tQ  =  X^/Uq.  The  conditions  required  for 
similar  solutions  of  the  conservation  equations  are  thus  renumbered  and  summarized  in 
Table  I. 

Table  I.  Conditions  required  for  similar  solutions  of  the  conservation  equations.* 


_ Similarity  condition 


(1) 

1- 

o 

0 

o 

1- 

o 

Mq 

(2) 

11 

o 

(3) 

®w0  ""  ®0 

(4) 

•o  Uo 

(5) 

(6) 


U 


=  1 


_ Conditions  imposed  on  two  flows _ 

The  ratio  of  "change  in  mass/ mass"  must  be  identical  in  the 
two  flows  during  the  observation  intervals. 


The  variables  describing  the  incoming  material  must  use  the 
same  scale  factors  as  the  main -flow  variables. 


The  time  must  scale  as  the  ratio  of  "distance/ velocity" 
scaling  factors.  This  condition  may  be  thought  of  as  linking 
the  geometric  and  dynamic  scaling. 


The  kinetic  and  internal  energies  must  scale  the  same.  (For 
two  ideal  gases,  the  Mach  numbers  must  be  identical.) 

For  pipe  flows,  the  usual  conditions  of  Reynolds -number 
similarity  enter  through  the  wall  shear  stresses. 


Pq^O  Pq  ^  . 
®o^o”®oP 


The  ratios  of  "pressure/ energy  per  unit  volume"  must  be  the 
same  for  both  flows.  This  condition  imposes  a  restriction  on 
the  equation  of  state  of  the  gases.  (For  two  ideal  gases,  the 
T*s  must  be  identical.) 


(8) 


^0^0 


=  1 


The  ratio  of  "heat  added  (or  lost)/internal  energy"  must  be 
identical  for  the  two  flows  during  the  observation  intervals. 


❖ 


Derived  through  the  use  of  auxiliary  condition  Mq  =  PqVq. 
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SECTION  II: 

AUXILIARY  EQUATIONS  USED  BY  PUFL 


4* 


The  eight  conditions  given  above  are  very  general.  They  make  no  assumptions 
about  the  auxiliary  equations  for  mass  flux,  energy  lost  (or  added),  equation  of  state, 
wall  shear  stresses,  or  pipe  geometry.  The  auxiliary  equations  used  by  the  PUFL 
calculation  must  now  be  examined  in  order  to  relate  the  derived  similarity  conditions  to 
more  specific  physical  requirements. 


Equation  of  state 


First  consider  two  ideal  gases:  p  =  (7  -  l)ep  and  p'  =  (7  -  l)e'p  . 


Then, 

-  EL  -  (7  -  Deo  -  (7  -  1)  -  . 

PO  ■  ^  ■  (7'  -  De'  p’  (7’  -  1)  ®0P0- 

Similarity  condition  (7)  requires  that  Pq/^qPq  ~  hence,  for  the  two  flows  to  be 
similar,  it  is  required  that  7=7. 

For  calculational  purposes,  it  is  usually  convenient  to  express  the  equation  of 
state  as  a  polynomial.  The  following  example  is  used  here  to  indicate  that  condition  (7) 
is  really  a  little  more  general  than  the  frequently  stated  similarity  condition  which  re¬ 
quires  that  7,  the  ratio  of  the  specific  heats  of  the  two  gases,  be  constant  and  equal. 

Suppose  that 

p  =  +  a2p  +  a^e  +  a^ep  +  .  .  . 


and  that 


j. 

0 


'  =  -E-  =  +  +  -= —  +  .. 


P  ■  Po  '  Po  '  Po  '  P 


Po 


But  similarity  condition  (7)  requires  that  Pq  =  e^pQ  =  ep/e'p  .  Hence, 
a. 


®oPo  ®o  Pq 

and  similarity  can  exist  for  gases  with  polynomial  equations  of  state,  providing  that  the 
polynomial  coefficients  are  related  as  shown  above;  namely. 


a,  = 


a,. 

^0 


1  e^pQ  -  ^0 

Note  that  ideal  gases  are  a  special  case  of  the  polynomial  with  a^^  -  8-2  ~  Sg  ”  Oj  and  with 


^4  =  7. 
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Wall  Shear  Stress 

3 

The  frictional  stress  (T  )  exerted  by  the  wall  is  shown  by  Schlichting  to  be 
o  ^ 

proportional  to  pu  .  The  proportionality  may  be  replaced  by  an  equality,  * 

=  Cj  (Re,  e/D)  pu^,  by  using  the  dimensionless  coefficient  of  friction  (C^,  which  is 
a  function  of  the  Reynolds  number  (Re)  and  the  ratio  of  surface  roughness  to  diameter 
(e/D).  For  flows  at  low  Reynolds  numbers,  is  nearly  independent  of  surface 
roughness:  C^-  Cj(Re).  For  turbulent  flows  at  high  Reynolds  numbers,  is  virtually 
independent  of  the  Reynolds  number  but  depends  strongly  on  the  surface  roughness: 

Cf  —  C^(e/D).  Because  ablation  requires  high-energy  flows  that  are  usually  well  into 
the  turbulent  range,  for  such  flows  the  ratio  e/D  is  of  prime  importance. 

If  CjQ  =  Cj/Cj,  then  =  C^qPqUq,  and  similarity  condition  (6)  may  be  written 
as 

^wO®o^o  _  ^foPo^o^oV  -  ^fo^o^o^o  =  ^fO^O^O  ^  , 

UqMo  -  UoPqVo  -  Vq  "  Vp  -  • 

Consider,  for  instance,  two  straight  pipes.  In  this  case,  Sq/Vq  =  Then, 

condition  (6)  requires  that  =  Rq/Xq.  Suppose  that  the  axial  distances  scale  by  a 
factor  of  10  (Xq  =  10),  and  the  radial  distances  by  a  factor  of  20  (Rq  =  20).  Then, 
similarity  requires  that  C^q  =  “  20/10  and  that  Thus,  the  friction 

coefficient  of  the  unprimed,  or  larger -radius,  pipe  must  be  twice  as  large  as  the  primed, 
or  smaller-radius,  pipe  for  the  flows  in  these  two  pipes  to  be  similar. 


Ablation  Equations 


For  dynamic  gas  flows,  the  turbulent  convective  heat  flux  may  be  written  as 
CrrPuh,  where  Crr  is  the  dimensionless  coefficient  of  heat  transfer  and  h  is  the  specific 
total  enthalpy.  The  rate  that  energy  is  lost  by  turbulent  convection  from  the  gas  flow 
to  walls  having  a  surface  area  S  is  then  C^^puhS  =  -H. 

Similarity  condition  (8)  may  be  written  as 


^0^0  CpjpuhS 

V^^C^P'u'h'^ 


®o^o 


^0 


2  t 

Using  similarity  conditions  (5)  and  (7),  it  can  be  shown  that  h  =  e  +  p/p  +  u  /2  =  h  Op. 
Also,  using  condition  (4)  and  Mp  =  PpVp,  and  introducing  the  expression  ^HO  =  Vh- 
similarity  condition  (8)  may  be  written  as 
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♦ 


,  % 

'Ho'  V. 


=  1. 


The  energy  required  to  ablate  a  unit  mass  of  wall  material  may  be  written  as 
E  +  rjh,  where  is  the  specific  energy  reqviired  to  vaporize  the  wall  material  and  n 
is  a  turbulent -transpiration  coefficient.'^  The  mass  flux  may  then  be  written  as 


Cjjpuh 

""  "  E^+nh 

where  heat  is  transferred  to  the  wall  by  turbulent  convection.  When  this  expression 
is  usedj  similarity  condition  (1)  is  written  as 


"mT"" ^ 


Vo 

M„ 


=  1. 


For  two  pipes  made  of  the  same  material,  the  assumptions  "  E^/E^  1  and 

t)q  =  n/n'  =  1  may  not  be  too  gross  to  make.  Then,  because  h  =  h'  from  conditions 
(5)  and  (7),  and  using  condition  (4),  similarity  condition  (1)  reduces  to 

^0^0 

Vq 

with  the  auxiliary  conditions  that  E^q  =  1  and  tJq  =  1.  Hence,  conditions  (8)  and  (1) 
reduce  to  the  same  requirement. 

For  turbulent  flows,  Reynolds'  analogy,  which  assumes  that  the  same  mechanism 
causes  the  exchange  of  momentum  as  well  as  of  heat,  is  frequently  used.  Reynolds' 
analogy  may  be  expressed  as  Cjj  =  0^/2,  and  in  this  case  as  ^HO  ^  ^f0‘  Hence,  if 
Reynolds'  analogy  is  assumed  to  hold,  similarity  condition  (6)  becomes  superfluous 
providing  that  condition  (1)  is  satisfied. 

When  one  is  trying  to  experimentally  scale  a  flow  during  which  ablated  wall  ma¬ 
terial  adds  significantly  to  the  original  mass  of  the  flow,  similarity  conditions  (2)  and 
(3)  become  increasingly  important.  In  discussing  similarity  condition  (1)  above,  it  was 
found  that  the  specific  heats  of  ablation  for  the  two  flows  had  to  be  equal 
^Ey  +  T)h  =  E^  +  Ti'  h'j,  which  implies  the  same  pipe  material.  Hence,  e^Q  «  1 
probably  exists  under  most  such  conditions.  Also,  e^Q=  1  implies  Bq  =  1,  and  condition 
(5)  then  implies  that  Uq  =  1.  Thus,  if  significant  mass  addition  occurs  by  natural 
ablation  of  wall  material  during  a  flow,  a  scaled  flow  must  have  the  same  velocity  and 
specific  internal  energy  (temperature)  as  the  original  flow. 

If  the  added  mass  is  not  a  significant  fraction  of  the  original  mass  during  the  flow 
period,  then  conditions  (2)  and  (3)  lose  importance.  In  this  case,  only  the  requirement 
that  Uq  =  eQ  is  important,  and  the  possibility  exists  of  using  both  conditions  at  a  lower 
velocity  and  temperature  flow  as  a  scale  model. 
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A  number  of  possibilities  exist  for  specifying  the  kinetic  and  internal  energies  of 

the  ablated  mass.  The  important  requirement  that  the  kinetic  and  internal  energy 

partition  for  the  entering  material  be  the  same  in  both  flows  must  hold;  i.  e., 

fi  =  U^/2e  and  =  1.  This  requirement  arises  from  similarity  conditions  (2),  (3), 
w'  w  '■^0 

and  (5). 

A  fairly  general  method  presently  used  by  PUFL  to  determine  e^  and  is  as 
follows.  The  rate  at  which  energy  enters  the  flow  is  some  fraction  (a)  of  the  rate  at 
which  energy  is  being  lost: 

e  +  — 

_  w  2 

“  E  +  ph  • 

V  ' 


The  fraction  a  may  be  determined  by  additional  equations  or  may  be  treated  as  a  con¬ 
stant  for  a  first  approximation.  For  o’  =  1,  there  is  no  total  energy  loss.  Using 
EyO  '  ^0  "  ^  conditions  (2)  and  (3),  similarity  requires  ofQ  =  e^Q,  which, 

as  discussed  earlier,  frequently  is  physically  required  to  be  unity. 

For  the  auxiliary  equations  used  by  PUFL  to  simulate  the  ablation  process,  the 
eight  previously  described  general  similarity  conditions  may  now  be  restated: 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


^HO^O^O 


UwO  =  Uo 


®w0 


"  ®v0  ’^0  ^  ^ 


or,  with  the  PUFL  equations, 
it  is  sufficient  that 


“o  ®w0 
^0  =  1 


u; 


^fO^O^O 

Vn 


=  1. 


If  Reynolds'  analogy  holds,  this  is  satisfied  by  condition  (1). 


Pq  =  ^0^0'  equation  of  state. 

Presents  no  new  requirements  if  the  requirements  of  conditions  (1)  are 
satisfied. 
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SECTION  III: 

CALCULATIONAL  VERIFICATION  OF  THE  SIMILAR  CONDITIONS 


In  order  to  test  the  previously  derived  conditions  for  similarity,  some  PUFL 
problems  are  now  considered.  First,  two  finite,  high-energy  slugs  of  gases  are 
considered  to  be  moving  into  lower-energy  gases  in  two  straight -walled  pipes  with 
frictional  effects,  energy  losses,  and  addition  of  mass.  The  initial  conditions  for  an 
unprimed  and  a  primed  flow  are  shown  in  Fig.  1. 


From  the  initial  conditions,  the  following  scale  factors  are  obtained: 

Pq  =  1.0  X  10^,  Pq  =  1.0  X  10^,  e^  =  1,  Ug  =  1,  Xg  =  10,  and  Rg  =  5.  Hence,  similarity 

conditions  (5)  and  (7)  are  satisfied.  The  two  pipes  are  straight  for  their  entire  length, 

hence  Sg/Vg  =  1/Rq-  Condition  (1)  then  requires  that  Cj^g  Xg/Rg  =  1  or  that 

C  =  R^/X„  =  1/2,  hence  C,,  =  C'  /2.  For  the  two  problems  considered  here, 

HO  0^  0  '  tL  hi 

Crr  =  0.002  and  =  0.004  have  been  picked.  In  addition,  E  =  E'  =  4.0  X  10  and 
H  hi  V  V 

Y]  =  T]'  =0.1  have  been  chosen.  It  is  assumed  that  one -half  of  the  energy  lost  from  the 
flow  is  returned  as  ablated  mass  (a  -  a  =  0.5),  and  that  ^  =  "  U^/ 2e^  to 


nal  energy  —  ergs/g 
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satisfy  conditions  (2)  and  (3).  Condition  (4)  is  satisfied  if  the  time  is  scaled  by  the 
factor  10  =  tQ  =  Xq/Uq.  Reynolds’  analogy  is  assumed  to  hold  (i.  e.,  C^jq  =  C^q),  and 
condition  (6)  is  satisfied  by  =  2C^  =  0.004  and  C|  =  2C^  =  0.008. 

The  results  of  the  puFL  calculations  for  these  unprimed  and  primed  flows  are 
shown  at  times  of  0.8  msec  and  0.08  msec  in  Figs.  2  and  3  respectively.  These 
figures  show  that  the  solutions  are  similar. 

The  next  two  problems  deal  with  pipes  having  non -constant  radii  and  converging 
nozzles.  The  surface  area  and  volume  for  the  frustums  of  cones  (i.  e.,  cylindrical 
nozzles)  are 


S  =  7r(r^  +  r^)  -  ^2^^ 


V  =  (r^  +  +  ^2)- 


Hence,  in  order  to  maintain  a  constant  ratio  of  Sq/Vq  =  1/Rq  for  nozzle  problems,  the 
axial  length  and  radii  must  scale  the  same  (i.  e.,  Xq  =  Rq). 


Unprimed  flow 


r  =  100/\/^ 


0  cm 


r  =  SO/nTF 


Primed  flow  (geometrically  scaled  by  a  factor  of  2) 


Fig.  4,  Initial  conditions  for  unprimed  and  primed  nozzle  flows  (flow  parameters  are 
the  same  as  those  in  Fig.  1). 
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The  two  flows  considered  next  are  again  finite,  high-energy  slugs  of  gases 
moving  into  lower -energy  gases,  just  as  in  the  last  problem.  The  geometry  of  the 
unprimed  and  primed  flows  are  shown  in  Fig.  4.  This  geometry  imposes  the  scale 
factors  Rq  =  Xq  =  2.  From  the  initial  conditions,  Uq  =  1;  hence,  t^  =  Xq/Uq  =  2  =  t/t' . 

The  pipe  walls  are  allowed  to  radially  expand  a  distance  equal  to  their  original  radius. 

The  expansion  time  is  5  X  10~^  sec  for  the  unprimed  flow  and  2.5  X  10  ^  sec  for  the  primed 
flow.  Under  these  conditions,  at  any  instant  the  surface -to -volume  ratios  of  the  ex¬ 
panding  nozzles  remain  constant:  Sq/Vq  =  I/Rq  =  1/2.  These  problems  use 
Cj  =  0.008,  cj  =  0.008,  Cjj  =  0.004,  =  0.004,  along  with  the  same  ablation  conditions 

as  the  previous  problems.  Hence,  the  remaining  conditions  for  similarity  are  satisfied. 

The  results  of  the  PUFL  calculations  for  these  two  flows  are  shown  in  Fig.  5  (p.l6) 
for  the  unprimed  flow  at  a  time  of  0.6  msec  and  in  Fig.  6  (p.l7)  for  the  primed  flow  at  a  time 
of  0.3  msec.  These  results  are  shown  right  after  the  shock  starts  to  reflect  from  the 
rigid  right-hand  wall.  The  results  are  seen  to  be  similar,  indicating  the  validity  of 
the  derived  conditions. 


SUMMARY 


For  turbulent  flows,  the  familiar  Reynolds -number  requirement  for  similarity 
may  be  reduced  to  requirements  on  surface  roughness  and  surface -to-volume  ratios. 

For  heat  sources,  similarity  requires  that  during  the  scaled  observation  intervals, 
the  ratio  of  "heat  added/ internal  energy"  be  identical.  Similarity  for  mass  sources  re¬ 
quires  that  the  ratio  of  "change  in  mass/ mass"  be  identical  in  two  flows  during  the 
scaled  observation  intervals.  Also,  variables  describing  entering  mass  must  use  the 
same  scale  factors  as  the  main -flow  variables.  For  gas  flows  where  ablation  is  not 
important,  equal  Mach  numbers  are  required  for  similarity.  However,  if  significant 
mass  entrainment  occurs  during  the  observed  time  intervals  of  the  two  flows,  and  if 
the  entrained  mass  has  identical  properties  when  it  enters  the  two  flows,  then 
similarity  requires  that  the  two  observed  main  flows  have  the  same  velocity  and  temper¬ 
ature.  Geometrically  similar  turbulent  flows,  in  which  heat  transfer  and  friction  are 
related  by  Reynolds'  analogy  and  in  which  equivalent  assumptions  are  made  about  the 
entering  mass,  are  shown  to  readily  satisfy  the  derived  similarity  conditions. 
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APPENDIX  A 


The  approach  to  similarity  used  in  this  report  considers  two  arbitrary  physical 
flows,  the  primed  and  the  unprimed.  At  some  instant  in  time,  the  two  flows  are  exam¬ 
ined  and  ratios  of  all  of  the  pertinent  variables  are  obtained.  The  requirements  imposed  on 
these  ratios  are  then  determined.  These  requirements  must  be  satisfied  in  order  for  the 
differential  conservation  equations  to  be  identical  in  both  flows. 

This  approach  differs  slightly  from  the  approach  traditionally  used.  It  has  been 
chosen,  however,  because  of  its  slightly  closer  relationship  to  the  frequently  encountered 
physical  situation  of  having  two  experiments  and  asking,  "What  relationships  must  exist 
between  the  parameters  in  the  two  experiments  at  this  instant  in  order  that  one  may 
represent  a  scaled  model  of  the  other  as  time  increases?"  It  is  shown  in  the  following 
discussion  that  the  approach  used  in  this  report  and  the  one  traditionally  encountered  in 
similarity  analysis  give  equivalent  results. 

For  the  more  traditional  approach,  dimensionless  quantities  are  obtained  from 
the  flow  parameters  by  selecting  certain  suitable  and  characteristic  magnitudes  for  the 
flow  being  considered.  For  example,  consider  a  pipe  flow  described  by  a  set  of  un¬ 
primed  variables.  A  characteristic  length  of  the  flow  (Xj^)  is  s^et  equal  to  the  radius, 
and  a  reference  time  (tj^)  is  frequently  obtained  from  a  reference  velocity  (Uj^)  as 

^  reference  density  (p^),  mass  flux  (Mj^),  surface  area  (Sj^),  and  volume 
(Vo)  are  also  selected.  A  set  of  dimensionless  variables,  denoted  by  asterisks, is  then 
obtained  from  the  physical  flow  variables  and  from  the  reference  magnitudes: 


X 


X. 


R 


m 

•  >  ^ 


and  p*  =  . 

pR 


Derivatives  in  the  dimensionless  system  are  related  to  the  derivatives  for  the  flow 
variables  as  follows: 


a 

=  Q 

8X 

1 

a 

ax 

9X 

ax 

■^r‘ 

❖ 

ax 

D 

D 

Dt"  1 

,  D 

Dt 

"  Dt  ET  ' 

Dt 

‘'R 

Dt 

and 
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The  continuity  equation  for  the  flow  variables  is 

Dp  _  jO  .  ^  +  Sm 
Dt  ■  ■ V  dt  V  • 

Using  the  relationships  defined  above,  the  continuity  equation  becomes 


_  * 

‘h“* 


^  • 

p  Pj^dV  S  Sj^m  noj^ 

-I- 

V  Vj^tj^dt  V  Vj, 


Hence, 


❖ 

»i< 

Dt 


dV 

,  * 
dt 


^R^R 


S  m 
>[< 
V 


is  the  dimensionless  form  of  the  continuity  equation  for  the  flow. 


Now  consider  another  flow  denoted  by  a  set  of  primed  variables.  A  dimensionless 
set  of  variables  are  obtained  from  the  primed  flow  variables  by  using  a  primed  set  of 
reference  values.  Following  the  same  procedure  as  that  used  for  the  unprimed  flow, 
the  dimensionless  continuity  equation  for  the  primed  flow  is 


Dt 


I  5JC 


dV 


dt 


I 


f  .  I  I 

^r”^r*r 


S  m 


V 


In  order  for  the  primed  and  unprimed  flows  to  be  dynamically  similar,  the  primed 
and  unprimed  dimensionless  continuity  equations  must  be  identical.  Hence,  dynamic 
similarity  requires  that 


The  reference  values  may  be  chosen  arbitrarily.  Hence,  no  loss  of  generality 
results  if  the  reference  values  are  chosen  so  that  all  of  the  primed  and  unprimed 
dimensionless  variables  are  equal  to  unity.  Then,  for  instance, 

S*  =  s 

«  Sj, 


and 


S  _  ®R 


^R 


Hence,  the  ratio  of  the  variables  in  the  two  flows  can  be  expressed  in  terms  of 
the  ratio  of  the  reference  values.  This  report  gives  the  ratio  of  reference  values  as 

.1  t 

Sp^/S^  =  Sq  =  S/S  .  Thus,  the  notation  used  in  the  report  may  be  expressed  as 
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which  is  the  same  as  that  obtained  in  the  report. 

Generally,  the  traditional  approach  of  determining  whether  two  flows  are  similar 
involves  comparing  groups  of  reference  values  like  "that  are  derived  from 

continuity.  When  all  such  groups  (from  the  complete  set  of  equations)  are  numerically 
equal  for  two  flows,  the  flows  are  called  similar.  The  flows  then  behave  in  a  similar 
manner  as  time  advances. 


Although  the  reference  values  have  dimensions,  the  referred-to  groups  of  these 
reference  values  are  dimensionless.  Their  nondimens ional  character  is  due  to  the  fact 
that  they  are  coefficients  of  terms  in  the  dimensionless  set  of  equations.  For  instance, 
the  Reynolds  number  that  is  frequently  used  in  similarity  work  is  just  such  a  dimension¬ 
less  group  obtained  from  certain  forms  of  the  momentum  equation. 

For  the  PUFL  equations,  it  is  not  intuitively  obvious  what  values  should  be  chosen 
for  reference  or  as  characteristic  values  for  some  variables  (e.g.,  the  mass  flux). 

Hence,  it  is  legitimate  to  examine  both  flows  at  some  instant  and  to  choose  the  reference 
values  so  that  the  dimensionless  variables  have,  at  that  instant,  a  numerical  value  of 
unity.  Then  the  ratio  of  reference  values  becomes  simply  the  ratio  of  flow  variables. 
This  is  the  essence  of  what  is  done  in  this  report. 
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APPENDIX  B 


The  momentum  and  energy  equations  given  in  this  report  are  examined  in  the 
following  discussion,  and  conditions  for  similar  solutions  to  the  equations  are  derived. 


Momentum 


The  momentum  equation  is 


ir  =  m 


Substituting  for  the  primed  variables  gives 


^^0  21L  -  1  fm'  M  S'  S  ni  u'  -  U  u' )  -  ^ 

tp  ■  Dt'  'm'MQ  ®0'^w0'^w  ^0  '  rn'M^X^  dx 

The  energy  equation  in  the  primed  system  thus  is 


,  ^O^'Po  9^  v 

'  "  ni'M^Xg  ■  8x'  m'  Mq 


Du'  ^0  ^O^O^wO  rh'S',/  m'S',,  _^0  ^0^0  v'  8p’ 

dT'U^"  Mq  ■  m'  w  UqMq  m'  Uq  ‘  XqMq  m'  9x' 

*0  '^wO^O  ^w^ 

■  Uq  •  Mq  •  m'  • 

From  the  above,  the  conditions  required  for  similar  solutions  are: 

^0^0®0^w0_, 

MqUo 

(3)  joyopp-  =  1 

UoXqMo 

^o’^wO^O  _  , 

- ^ 

*^0  0 

Since  condition  (2)  requires  that  ^Q^QtQ/MQ  =  1,  condition  (1)  imposes  the 
additional  requirement  that  U^q/Uq  =  1. 
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Energy 


The  energy  equation  is 


De^  mS  ,  mS 


Dt  m 


+  —  (e  .  e)  -  + 

m  '  w 


.  .  J2Y  , 

m  *  Dt  m 


m 


Substituting  for  the  primed  variables  gives 


12.  .  ^ 

*0 


^0^0  m'  S 


M. 


(“'Uo-“;u»o>  *o”'s'So  ,  ,  ,  , 

m'  '  2  m'  Mq  ^®w®w0  "  ®  0^ 


PoP'  I^„I^,SoS'|u'|  H'H„ 

U  A 


M^m 


T  • 


t„  dr 


m'  Mq 


0  m'  Mo  ‘ 


The  energy  equation  in  the  primed  system  thus  is 


De'  _  ^0  „2  m’s’  u' 

Dt'  "  Bq  •  Mq  ^0  ni'  •  2 


'"0 


M. 


UwO^O 


m'S'u'u 


w 


m 


»2 


^  2  m's' 

Gq  ■  Mq  ^0  m'  '  2  CqMq  wO 


m'S'u'u;^ 

I 

m 


_  ^  m’s'  e'  _  ^pPp^O  _  p'pv' 


epMp 


m 


®P^P^P  m'Dt' 


VwoSoUo  |u'|  _ 


^0^0 


m 


eoMom' 


The  conditions  for  similar  solutions  imposed  by  the  conservation  of  energy  thus 


are: 


(1) 


^p^p^p^p 

®p^o 


(2) 


^P^P^P^P^WP 

®P^P 


(3) 

(4) 


^P^pSp^wP 

®P^P 


=  1 


tpMpSpe^P  ^ 

"P^P 
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(5) 


Po^o 
iT  =  1 
®o  0 


.  >  VwO^O^O  _ 


*0^0 
®0  0 


Since  condition  (5)  requires  that  MQSQt^  =  M^,  conditions  (1)  through  (4)  may  be 
simplified  to: 


UqUwO 


_w£  =  1 
®o 
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